On the density of the winding number of planar 
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Abstract 

We obtain a formula for the density f{0,t) of the winding number 
of a planar Brownian motion Zt around the origin. From this formula 
we deduce an expansion for f{\og{^/i) 9 , t) in inverse powers of log y/i 
and (1 + 0^)^/^ which in particular yields the corrections of any order 
to Spitzer's asymptotic law ( [H]). We also obtain an expansion for 
f{6, t) in inverse powers of log ^/i, which yields precise asymptotics as 
t oo for a local limit theorem for the windings. 

Keywords: planar Brownian motion, winding number, transition density, Spitzer's 
law, local limit theorem, asymptotic expansions 

1 Introduction 

In his celebrated paper about 2-dimensional Brownian motion, F. Spitzer 
considered the transition probabilities of &t, its continuous winding number 
around the origin, and showed in particular that the limiting distribution of 
, is the Cauchy distribution: 

logv^ 

limPf i- < a) = - / -. (1.1 

Several proofs of this result have since then been given, see for instance 
page 43 of the book [lOj by L.G.C. Rogers and D. Williams and the discussion 
and references therein. Extensions to more general situations are known as 
well, for instance those considered by J. Pitman and M. Yor in [2], [5] and [U] 
in the frame of asymptotic laws for planar Brownian motion. 

More recently, V. Bentkus, G.Pap and M. Yor in ( [Ij and [6j) obtained 
by Fourier methods an expansion of the distribution function of G^, that 
yields detailed asymptotics, both in t and a, as they go to oo. 
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F. Delbaen, E. Kowalski and A. Nikeghbali proved in [3J (also using 
Fourier methods) the foUowing local limit theorem for Qf. if a < /3 G M, 

lim log\/tP(a < et < ;5) = (1.2) 

In the present work, we obtain a formula for the density f{0,t; p) of 0j 
when the initial condition Zq = p 7^ 0, p G M, that follows by integrating an 
expression for the joint density p{r,9,t;p) of {\Zt\ ,Gf) that appears in [2j 
(and follows as well from results in [12], see also the comments there for 
older references). We compute / in the next section. Let us mention that 
in Chapter V of [5], R. Mansuy and M. Yor discuss some representations for 
the distribution of the winding number of the Brownian lace of lengh t. 

The expression for the density / is given in terms of an integral including 
a couple of fractions. In the last section, we expand these fractions in inverse 
powers of log \/t and obtain two asymptotic expansions, after term by term 
integration. The first one is an expansion for f (log \/t 6 , t; p) very close to 
that given in [5]. The difference is that the dependence of the coefficients on 
6 is given more explicitly, and that we obtain also accurate estimates on the 
behaviour in t of those coefficients. A more precise comparison is discussed 
in Remark 4 after the statement of the corresponding result (Theorem 13. ip . 
The second expansion is for f{6, t; 1) in inverse powers of log \/t, that yields 
in particular the corrections of any order to (jl.2p . The result is stated as 
Theorem [321 



2 A formula for the density 

We consider a planar Brownian motion Zt starting at zq 7^ 0, 0t a continuous 
determination of its argument and Rt = \Zt\. Let us call p the joint density 
of Gt and Rt, 

p{r,9,t; p,a)rdrd0 = F{Rt € dr,@t G d9\Ro = p,Qo = a) 

Since p is clearly a function of 9 — a, it suffices to consider a = 0, what 
we do in the sequel. The following expression for p was deduced in [2] (see 
also [12]), 

p{r, 9, t;p) = — exp {-^-^^) T cos{i^9) I A du (2.1) 
Integration in r of the above formula yields an expression for the density of 

e*, f{9,t-p) d9 = Pieted9\R^ = p). 
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Denote by 1a the indicator function of the set A. We have then: 
Proposition 2.1. The density f of Qt is given by 

f{0,i;P)= 7= cos e I(_^/2,7r/2) {^) 

\/ Zirt 

oo 

+ f e-^"°"^'^sinh(a;/2) (2.2) 

27r v27rt j o 

. w/2 u/2 . 

^ H^^/2)2 +(0 + 7r/2)2 + (w/2)2 +(9-tt/2)'^'' 

Before proving the previous proposition, let us show another formula 
for the density /, that will prove more suitable to obtaining asymptotic 
expansions (as f — )• oo). It follows directly from the change of variables 

2 

^ coshcj — )• z in the integral above. 

Corollary 2.2. The following formula for f holds: 
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-^(1-cos 261) 



p e it ' 

f{0,t; P) = 7= COs6II(_^/2,^/2)(6') 

\ lilt 

+ 7r-^ / X (2.3) 

iarcosh(^) iarcosh(^) 
^ (i arcosh(^))2 + (0 + ^/2)2 (i arcosh(^))2 + (0 - 7r/2)2 ^ 

Proof of Proposition \2.1\ Integrate by parts once, recall that 2/^(x) = 
Iy-i{x) + Iyj^i{x) and use formula 6.618-4 of [5] to compute 

e~2tl^[ — )rdr = p I e~2t/^( — )dr = 

e -[/.-i(^) + W(^)]dr = ^e..[/^(^) + /^(^)] 

With the aid of this formula we can integrate (j2.ip . after interchanging the 
order of the integrals, obtaining 



oo 



f{e,t;p) = I p{r,e,t;p)rdr 


2 

P j cos(^0)[/^(^) + /^(^)]d^ (2.4) 



2\/27rt 
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Recall next the integral representation for the Bessel function 1^ (formula 
8.431-5 of m): 

71" Jo Jo 
Substitution of this last formula in the integral in (I2.4p yields 



cos(z.e)[l^(^)+/^(^)] dv 



1 roo i-K 2 

-/ cos{u9) e'^^"'"^ [cos{{u-l)uj/2)+cos{{u + l)uj/2)]dujdu 
^ Jo Jo 

1 rco poo 2 

/ cosH)sin(^(i/- l)/2) / e-iF™«hM-('^-iW2^a;dz^ 

TT Jo Jo 

fOO />oo 2 



■y poo pc 

/ cos(i/6l)sin(7r(i/ + l)/2) / 

71" Jo Jo 



^1 + ^2 + ^3 (2.5) 



To compute the first integral ^i, use that cos((z^ — l)a;/2)+cos((z^+l)a;/2) = 
2cos{i'uj/2) cos{uj/2), and recall that the integrand is an even function of oj 
and 1^ to write: 



1 f°° f2 
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S^ = — I cos(z^0) / " 2e« ^'^^.^coso; cos{v uj) dujdv 
2vr J_oo J-f 

= 2e«^°"2epog^ (2.6) 

The last identity follows from the Fourier inversion theorem, which applies 
since the inner integral is an function. This in turn can be seen from 
the fact that both and 1 1^^/(1/) | are uniformly bounded. The last state- 

ment can be checked integrating by parts twice, using that 2 e ^ cos uj 
vanishes at it^. 

To compute 52, change the order of the integrals, recall that 
sin(7r(zv — l)/2) = — cos(z^7r/2) and integrate 



roo 

\ COS (vQ) cos (z^7r/2)e ~ dv 
Jo 



2^(y/2)2 + (^ + vr/2)2 (y/2)2 + (0 _ vr/2)^ 
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The computation of 6*3 is similar, and we obtain finally 



1 /'oo 2 

^2 + ^3 = - / e-^'"'^'^ sinh(a;/2) 
71" Jo 



/ uj/2 uj/2 \ 

V (a;/2)2 + {e + ^/2)2 + + {6- 7t/2)^ ) 

(2.7) 

From ([23]) . ([M]) and this last equation yield ([22]). □ 



Despite formulae (|2.2p and (j2.3p may appear complicated, we can see 
that they give information on the behaviour of Qf In particular, it is not 
hard to deduce for instance detailed asymptotics as i — cx). We state and 
prove two precise results in the following section. 

3 Two asymptotic results 

Theorem 3.1. The density f{6,t;p) admits the following expansion: for 
any natural number N > 0, t > 2 and 5 G (0, 

logiV-t)f{9log{V-t),t-p) = (-ir-4ng)Cn(t;p) ^ 

2vrA^o 2''-^{l + 9^) — {log{Vt)) 

(3.1) 

The coefficients satisfy 

Y^keven ("fc ^) (" 1) ^ 

An{e) = ^^ , Cn{t;p)=Cn + On{t-') (3.2) 

(1 + 02) — 

where the Cn are constants given by 

(\ /"OO — z 

^{log{8z/p')r-'dz. (3.3) 

fc <n 

The remainder terms Rn satisfy that for some positive constants k^, 

\Rn\ < — — (3.4) 

(log(^/t))^+i(l + ^')^+' 



Let us make some observations on the previous result. 
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Remarks 



(1) The dependence on the initial condition p is made exphcit in the co- 
efficients Cni but the constants /cat, and On{t~^) are not uniform in 
p. We could have taken p = 1 as usual, but preferred to keep it ex- 
plicit, as the dependence may be traced back in order to investigate 
the behaviour in p also. 

(2) For N even, ~ , ^2-11/2 as — )• 00, while for odd, Aj^ ~ 



(1+6*2)1/2 

1 as — )• 00. Then the order in 6 of the coefficient — "^^^^n+i of 



{\ogVtr 



(1+^2) 

in the expansion is the same as that of Rn for N even. 



(3) As Co = \/vr and ^0(6*) = ^j^_|_g2-)i/2 > the first term yields precisely 
Spitzer's result (jl.ip . The c„ can be expressed in terms of derivatives 
of the Gamma function T, using the formula 

^(log(.))"d. = |lr(«)|,^, (3^5) 

(see for instance formula 4.358.5 on page 578 of [4]). 

(4) As mentioned in the introduction, V. Bentkus, G.Pap and M. Yor in 
( ^ and ^j) obtained an expansion for / (in the case p = 1), by 
different methods. The expression for the coefficients in 6 in (j3.ip is 
slightly more explicit here than the one they obtain. It is not difficult 
to see that the first three terms coincide with those they compute. 
We also provide accurate estimates on the rate of convergence to of 
Cn — Cn as t ^ 00. 

As the density / of Gt is a function of we consider p = 1 in the next 
result, and omit the reference to the initial condition, keeping the notation 
f{0,t). 

Theorem 3.2. For any natural number N > and t > 2, the following 
expansion holds: 

logiVt) f{e,t) = ^jz TTT^ + ^^-^^ 

27rVvr^ (log(Vi)) 
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The Qnid) are polynomials in 6 given by 

9n{e) = 

keven V 

fe <n 

(3.7) 

and the remainder Ejsj satisfies 

I p I 'S'jV+2(l^l) o^ 

(log(Vi)) 

Sisf-\-2 above is a polynomial of degree N + 2, whose coefficients may depend 
on N. 

Remark It is easy to compute go = 2-v/vr, and then the local limit theorem 
(jl.2p follows at once from Theorem 13. 2i Moreover, corrections of any order 
to (II. 2p in inverse powers of log ^/t follow from ()3.6p . Simple integration of 
the polynomials provide the coefficients, let us say G„, for the expansion 
of log-v/?P(a < 0* < /3) = X]n>o (iJ^^^Ti)^- instance, computing with 
the aid of (1331) and (ET 



gi = -2 ^(log 8 + log z) dz = 2V^(7 + log -) , 
Jo ^ 

where 7 is the Euler constant, we have 

log v/JP(o < e, < ,3) = ^ + -i^ (^-a)h + logl) ^^ , 

vr logvi 



The proofs of both theorems are given after stating and proving the next 
three lemmas. 

Let us shorthand 
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Lemma 3.3. The density f satisfies: 

1 /-oo -(^+p2/4t) 

logiVt) f{e log(^/t), t; p) = M{e, t;p) + 



27r0r Jp2/u y/z + p^/At 

( l + bx 1 + hx \ ri a ^(\\ 

V(l + 6x)2 + (6' + fx)2 + (l + 6x)2 + (0-fx)2; ' ^ ^ ^ 

T/ie ierm A/" is negligible as t — t- oo, uniformly in 9 in the following 
sense: for any given K £ N and 6 < 1/2, there exists a positive constant C, 
independent of p and t such that 

sup{l + e^)^ \M{e,t-p)\ < Cpt-^ for anyt>2 (3.11) 



Proof. Consider the expression for log{\/t) f {6 log{\/t) , t; p) obtained from 
()2.3p . and denote by the term coming from the first term on the r.h.s.. 
It is clear that it satisfies (jS.lip . 

The other term comes simply from organizing the fractions in the inte- 
gral, (jS.gp and the formula arcosh y = log [y + \/ y'^ — 1) ■ □ 

In the next lemmas, we show how to obtain an explicit expansion in 
powers of x of the fractions in the last line in (|3.10p . Let us denote by 3f?(z) 
the real part of z G Z. 

Lemma 3.4. Consider, for b,c,d G M such that (? — Ad < 0, the rational 
function 

. . 1 + bx 

' 1 + cx + dx^ 

For any given natural N > 0, the following finite N -expansion with remain- 
der holds: 

N 



9{x) = ^anx" + qN{x), 



n=0 



C— 11" 2b — r , 

where a„ = K((l + i^==)(c + iv^id^r) n>-l 

(3.12) 



j^j^i daN-i + caN + xdaN . 



q^ix) = -x^^+^( — — ^ — ) (3.13) 
1 + cx + dx'^ 
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Proof. The sequence a„ has to satisfy the recurrence an+2 + ca„+i + da„ = 
for n > 0, with gq = 1 and ai = b — c. To solve it, let us introduce the 
function ^ 

n] 



F(x) = ^^ (3.14) 



n>0 



From the recurrence for the a-n's, it follows that F has to satisfy the dif- 
ferential equation F" + cF' + = 0, with boundary conditions F{0) = 
1, F'(0) = 6 — c, whose solution is 

F(x) = ^e™ + ^e*^^ (3.15) 



with A = i + and r = _£±ly^^ZiE, Expanding the right hand side 

of ()3.15p in its Taylor series, and equating coefficients with the expression 
([3ll]l . we get a„ = 25R(^r"), which is exactly ([3T2]l . 

The formula (|3.13|) for the remainder is proven easily by induction for 

> 1. Indeed, compute 

, , N 2^ d + c(b - c) + xd(b - c) . 

g{x) -ao-{b- c)x = -x^{ \ \ ^ 

1 + cx + ax'^ 

to prove it for = 1. Suppose it true for A^, and write next 

N 

9{x) - ^ flnx" - aN+ix^~^^ = qN{x) - aN+ix^~^^. 

n=0 

Use the recurrence relation for the a„ to conclude that the r.h.s is qn+i- 

A direct computation shows that the formula is also valid for = 0, 
with a_i given by taking n = — 1 in p.l2p . □ 



Lemma 3.5. The fractions on the second line of (j3.10p can be expanded in 
powers of x as follows: for any given integer N > and x < A, 

1 + bx 1 + bx 

+ 



(l + 6x)2 + (e + |x)2 (i + 6a;)2 + (^- |x)2 
^ (-l)"P„(b) cos((n + l)/3) 

> -T7+T X + Qn, 3.16) 

^0 2-1(1 + ^2) — 



where 



e^^ = ^±^, P^{b)= J2 Q(26)-'=(-l)lv^^ (3.17) 

k<n,k even 



(1 + 02)2 
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and Qn = jv — ^NibjO) with 

(1 + 02)^ + 1 

\£N{b, e)\<K (____ + I|e|>2o) (3.18) 

for some constant K >Q independent of x, b,N and 6. 
Proof. Let > given, and set 

l + hx 1 , l + hx 



(l + 6x)2 + (0 + |x)2 i + e^'-l + cx + dx'^ 



), (3.19) 



for c = and d = ^^^fj^ ■ Observe that c2 - 4(i = -^g^^, so the 

conditions of Lemma [3.41 are satisfied and we may apply formulae (|3.12p and 
p.l3p to obtain, after some algebraic manipulations, 

rT^^ = E«!."-" + ^!J'. with ,3.20) 

n=0 

^^^^^LK((l + ,e)"+i(26-.vr)") and (3.21) 

(1) ^ _ ^+1 (6^ + vr2/4)(aff_,+ xaff) + (25 + vrg)aff 

l + 02 + (26 + 7re)x + (62+7r2/4)x2 ^"^^ ^ 

Since the second fraction in (|3.16p is just the first one evaluated at —0, the 
expansion for the former is: 

' n=0 

where the following expression for On results from evaluating ()3.2ip at —0 
and conjugating the complex inside SR, 



(1 + 62;)2 + (0 - f x)2 - 1 + 02 ( E «n ^ ^" + 9iv' ) > (3-23) 



a. 



(2) 
n 



^^^5}?((l + i0)"+H2b + ^vr)") Biidqfie) = q^^\- 



2"(1 +02) 
From (ISTTP and (l3:2ni we have 

^ + ^ = 2n(i~+g2)n 3? ((1 + ((26 - ivr)" + (26 + .vr) 
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Then, (IXTBD follows from (i3J9]l . dOOjl and with Qtv = j^iq^'^ + 

g]y ). To complete the proof, we only need to show (j3.18p . From (|3.2ip . 



I < 



|(l + i0)n+i(26-OT)«| (452 ^ ^2) f 



n-l 5 

2 



2"(l + 02)n 2"(l+e2) 

and then from ()3.22p . there is a constant Ki > such that 

kP'|<."^'A-. ''"''+'')^"; , (3.24) 
where D is the denominator in (j3.22p . Recall that x < 4 to estimate 
D = l + 9^ + i2b + ne)x + (6^ + vrV4) = (1 + 9^) (^\^ + 

^ (1 + ^ ) ( l + g2 M<20 + I|^l>20 

> (1 + ^')(^^^^I|e|<20 + ^I|e|>2o), (3.25) 

SO, from (1221]), 

I (1), , ^ jv+i (462 + 7r2)f+i I|g|<2o 

2^-1(1 + 6I2)T Hl + ^a;)2 ' ' ^ 

The same estimate holds for |, and ()3.18p follows. □ 



Proof of Theorem \3.1\ 

The proof consists in substituting the fractions in (jS.lOp by the A'^- 
expansion (j3.16p and then integrate term by term, recalling that only b 
depends on the integration variable z. Let > be given, recall that 
X = jQg^^ ; and write then 

log(^/^)/(01og(^/^),^;p) = M{e,t;p) + 

1 E(-ir.- r ^^^»(')^--+ 



2^\/^ ^ 2—1(1 + 02)^ V/4t + pV4t 

1 /-oo e"(-'+''^/^*) 



/ , (3-26) 
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Call 

Cn= / -^===Pn{h)dz (3.27) 

1 Q-i^+p^/^i) 

R^=M{e,t;p) + —^ / , QNdz. (3.28) 

From (|3.17|) . cos ((n + 1)/?) = ^n(6'), so to prove Theorem 13.11 it is enough 
to show that Cn and Rn satisfy p.2p and (|3.4p respectively. 
Recall the definition of b in (13.91) , and set 



1 8z f'°° 
6o = hm 6 = - log(^) ^ c„ = lim C„ = / P„(6o) dz (3.29) 

as follows from ()3.17p and ()3.3p . The c„'s are clearly finite, and we have 

(-^n = Cn + / , : ^ dz — ^ dz 

Jpyu yjz + /)2/4t 7o 

(3.30) 



To estimate the integrals above, let us shorthand w = — yl — (fe)^) 
and note that < w < (? j^tz < 1/2 for z > p^/4t to obtain 

\h - 6o| = -- log (\-w) = - log (1 + ) <w<^^ (3.31) 

2 2 1 — w ozt 

Denote by A;„ positive constants depending only on n, that may change from 
line to line. From (|3.17p and the previous inequality, it follows that 

|i^„(6)-Pn(6o)| < A:„|6-6o|(|6or~' + l), (3.32) 
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Then, let < (5o < 1/2, and estimate with the aid of and ([331]), 

/■oo -(z+p2/4t) p-z 



2/4t y/z + p^/At 

/ |P„(6)-P„(6o)|ciz + 



T /'OO /'OO 2 2 

<^/ e-^z-3/2(|6ori + l)<iz+ / ^|P„(6o)|(^ + -^)dz 



1 I, /■! 



t J pi /At t Jpi/4t 

+ ^{ z-^-^^^dz + 1) < A:„t-(i/2)+5o (3.33) 

To obtain the last inequality, we have splitted the integrals according to 
z < 1 or z > 1, after observing that, for Sn a polynomial of degree n, 
|5'n(&o)| •z'^" ^ kn for < z < 1 . For z > 1, the integrals are clearly 
uniformly bounded. The estimation of the last term in ()3.30p is similar: 

fp^/At —z 

1/ ^Pnibo)dz\<knt-^2+^\ (3.34) 

Jo Vz 



which together with (|3.33p and (|3.30p concludes the proof of ()3.2p . 

It remains to show that Rm satisfies (j3.4p . Consider first the case \0\ > 
20. From ()3.1ip and (I3.18p . it is enough in this case to show that 

— - ^^-^^^ 



■■/At ^Jz + p^jAt 

which can be seen by the same procedure as that to estimate p.27p (just 
considering (46^ + vr^)^"^^ instead of Pn{b)). 

In the case \9\ < 20, let us split the integral in (|3.28p according to 
z G {p'^/At, />^/4\/t), or z > p^ jA^Tt- For the latter, observe that 



logi - logt 2 ^ ' 
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so from (|3.18p . again a procedure as that to estimate the integral in ()3.27p 
yields that this integral is bounded. For the case z G {p^ /^t , /9^/4\/t), use 
the expression for Qn coming from (j3.16p and take < (5i < | 



^ (-irP„(&) cos((n + l)/3) N 



n=0 



2"-i(l + 02)^ 



</ (- —)dz+kN Y\Pn{h)\dz 

< t^i+'^i (3.37) 



To obtain the last inequality above, start by estimating the first inte- 
gral in the previous line. Observe that 1 + bx = ^'''^°'^M4^*/p ) ^ ^^j^j change 
variables y = arcosh(^) to obtain: 

— =7 ; — -dz<logt / —— TTT^^ 

p2/it y/z + p'^/it{l + bx) Jp2 ^ arcosh ( ^ ) 

^ pbgt Sinhy ^ .arcosh(v7) ^-^^^^ ^ ^ p log t ^ ^ i 

i/t Jo yVcosh y yi/coshy ~ i/t 

(3.38) 

The last integral in ()3.37p is bounded by t"*"^*^^, as can be seen from es- 
timates similar to those used to estimate (|3.27p . what finishes the proof of 
Theorem EH □ 

Proof of Theorem \3.S\ Recall that p = 1 and x = 1/ log ^/i. From ()2.3p and 
(j3.9p . we have 

logiVt) f{e,t) = 

i^e-i(— )cos.I( ./,./,(.) + ' r 

+ /-I , , , (3-39) 



(l + 6x)2 + (0 + f)2x2 (1 + 63;)2 + (0 - f )2x2 

Let us denote by k positive constants that may change from line to line, 
and split the integral above according to 2; S (l/4t, 1/4-v/t) or z > l/4\/f. 
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Proceeding as in the estimation of ()3.38p . it follows that 



1 f ^ 

( l + hx 1 + bx \ , , _ 1 (o A(W 



Expanding the fractions in (j3.39p with the aid of Lemma 13.41 we have 
1 -t- /j-r 

(TWTI^Tf)^ =?/"-^"<^-^*) (3.41) 



where is evaluated at —6, and, for d = 6 + (6* + |) , 

0) = + i(| + 0))" (3.42) 



^+ _ _ AT+i dAM-i + 2Mjv + xdAN , . 

From (I3.42p . it is easy to see that 



\Ar,{h,±e)\<{h' + C-±efy. (3.44) 

For z > Yj^^ the denominators in are > (^)^, as follows from p.36p . 
Therefore 

|T±|<x^+i^62 + (|±0)2)^+\ 
and proceeding as in the estimation of p.35p . we obtain 

7= / (T^+T+)dz<kN 7 . (3-45) 

27rV^ Ji/^Vt ^/^TWt - '^(logVt)^+i ^ ^ 

On the other hand, from (j3.42p . it is easy to see that 

\An{b, ±9) - A„(6o, ±0)1 < |6 - 6o| K (|6or + 1^1" + l) (3-46) 

Next, set 

i-no -z 

— {{Anibo,e) +An{bo,-e))dz, 
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which from (j3.42p is precisely (j3.7p . Note that the gn s are polynomials in 
e of degree n. To finish the proof, from (1339]) . (I33Q|), ([3:iT]) and ([OS]) it is 
enough to show that for 6i <\, 

r {Anjb, 9) + -9)) e-(^+V4*) (yL„(6o, g) + A^jbo, -6)) e"^ 

-/ ^ ^ ^ '-^ dz < knil + W)t—^+^\ (3.47) 

Jo 

what follows from estimates similar to those in (13.33P and ()3.34p . with the 
aid of (iOIH and (IMjl . □ 
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